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Abstract 

Mixing in the system may provide a sensitive probe for new 
physics beyond the Standard Model (SM) but has so far eluded exper- 
imental observation. The SM predictions are typically small (< 10^"^) 
for the mixing parameters x, y which, in the absence of charge-parity 
(CP) symmetry violation, measure the mass (x = Am/T) and lifetime 
{y = Ar/2r) difference of the CP eigenstates in the D° system. The 
asymmetric _B-factory experiments BABAR and Belle open up the op- 
portunity of measuring x, y with unprecedented statistical precision 
and sample purities. Results from BABAR and Belle, and from CLEO 
are reviewed. 
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1 Introduction 



Mixing phenomena, i.e. the oscihation of a neutral meson into its corre- 
sponding anti-meson as a function of time, have been observed in the 
and systems pp. Up to now, oscillations in the system have eluded 
experimental detection. 

This is in keeping with the Standard Model (SM) prediction that mixing 
in the system should be substantially smaller than for the and 
[HE]. The SM prediction, as discussed below, is a consequence of - 
mixing being dominated by intermediate states with light d and s quarks, a 
feature unique to the system among the neutral mesons. 

Mixing phenomena are characterized by two dimensionless parameters 
X, y. When charge-parity (CP) symmetry holds, they correspond to the 
difference in mass (x = Am/T) and lifetime {y = Ar/2r) of the CP eigen- 
states in the neutral meson system. Calculations within the SM typically 
predict that |x|,|y| < 10~^ for the j^. Due to the smallness of the 
SM prediction, mixing in the system may provide a sensitive probe for 
physics beyond the SM. 

Independent of their discovery potential, precise measurements of x and 
y in the system are desirable because even SM predictions for x, y span 
several orders of magnitude The main challenge in SM calculations of 

- mixing is estimating the size of SU(3) flavor-symmetry breaking 
effects (see below) (Hj. 

In this article, measurements of the mixing parameters x, y in the 
system are reviewed. Results on CP violation are only discussed when they 
arise as integral part of the x, y measurements. 

The most recent results come from the asymmetric i?-factory experi- 
ments BABAR and Belle, in operation since 1999. BABAR, at the storage 
ring PEP-II of the Stanford Linear Accelerator Center (SLAC), USA, and 
Belle, at the KEK i?-factory in Tsukuba, Japan, have accumulated the 
largest currently available charm samples (integrated luminosity ~ 100 fb~^ 
by the end of 2002) and can be expected to continue data taking till the 
end of the decade. Hence, BABAR and Belle open up the opportunity of 
studying charm decays with unprecedented statistical precision. 

The first collider experiment to employ the methods for measuring x, y 
now in use at BABAR and Belle was the CLEO experiment at the Cornell 
Electron Storage Ring (CESR), USA, in operation until 1999. This review 
describes these methods and the results obtained by CLEO, BABAR and 
Belle. Wherever appropriate, comparisons are made with results from fixed- 
target experiments, notably from E791 and FOCUS, that were dedicated to 



3 



charm studies during the 1991/1992 (E791) and 1996/1997 (FOCUS) fixed- 
target runs at Fermilab. 



2 Mixing formalism 

Mesons are produced in strong and electromagnetic interactions as flavor 
eigenstates with a weh-defined quark content. The production Hamiltonian 
is of the form TCq = Ti strong + 'Hem- Its eigenstates are the flavor or inter- 
action eigenstates M^, M^. The CPT theorem requires that particle and 
anti-particle have the same mass and lifetime. Thus and correspond 
to identical eigenvalues uiq , 70 of Tio . 

Mesons are observed by way of their decays which are governed by the 
weak force. The eigenstates of the evolution Hamiltonian that is responsible 
for their decay, Ti. = TCq + Tiweak, are the mass or decay eigenstates Mi^2 
with the corresponding eigenvalues mi^2 and 71^2- They are the physically 
observable states that obey an exponential decay law with a slope given by 
the respective lifetime 71^2: 

|Mi,2(t)) = |Mi,2(0))e-'"^i-2*e-(^i-2/2)*. (1) 

As a consequence of the difference between production and evolution 
Hamiltonian, a sample of neutral mesons, produced, e.g., as a pure 
sample, evolves in time as a superposition of and states with time- 
dependent coefficients: 

\M{t))=CMo{t)\M')+oj^{t)\W). (2) 

The time evolution obeys Schrodinger's equation: 



(3) 

This equation represents two coupled differential equations. Its eigenstates 
are the physically observable eigenstates Mi^2- As a consequence of the non- 
zero off-diagonal elements in the Hamiltonian TC, the mass eigenstates Mi^2 
are obtained as linear superpositions of the flavor eigenstates: 

|Mi,2) =p|M°) ±g|MO). (4) 

In the absence of CP violation, i.e. for \q/p\ = 1, Mi^2 are CP eigenstates. 
For a more detailed discussion, see Ref. jS]. 
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2.1 The mixing parameters x, y 

Mixing phenomena can be characterized by two dimensionless parameters: 
y = (71 - 72)/(7i + 72) = Ar/2r, x = (mi - m2)/r = Am/r. (5) 

Mixing in a neutral meson system occurs when at least one of the two 
following possibilities applies: Between the decay eigenstates of the system 
there is a non-zero lifetime difference (y 7^ 0), or there is a non-zero mass 
difference (x 7^ 0). In the first case, flavor mixing is a consequence of the 
shorter- living decay eigenstate dying out. The remaining longer-living one is 
a linear combination of a and a component, so that even an initially 
pure sample shows some fraction of after some time. In the second 
case, flavor mixing is a consequence of a pure transition from to and 
vice versa. 

In the -D" system, according to the SM, CP violation is a negligible effect 
(see below); the decay eigenstates are, to a good approximation, also CP 
eigenstates. In the following, the eigenstate Di {D2) with mass mi (m2) and 
width 71 (72) is chosen as CP-even (CP-odd).^ In the K° and P° systems, 
where CP violation is experimentally observable [HI, the conventional choice 
for Ki and Bi [K2 and B2) is the heavier (lighter) eigenstate. 

In oscillation plots that show the fraction of flavor states in an ini- 
tially pure sample of as a function of time, the parameter x is related 
to the frequency of the oscillation, while y is related to the damping of its 
amplitude. 

The two parameters x, y reflect different mechanisms through which 
mixing can proceed in lowest order in the SM [7j. A can oscillate into 
its anti-particle by way of on-shell intermediate states that are accessible 
to both particle and anti-particle. These are long-range processes with am- 
plitudes oc —iy. On the other hand, a can also oscillate by way of 
off-shell intermediate states that can be represented by box-diagram loops 
(see Fig. These are short-range processes with amplitudes oc x. Ad- 
ditional contributions to the box diagrams in Fig. fright) from as-of-yet 
unobserved non-SM particles could result in a deviation of the measured 
value of X from the SM prediction. 

The time-integrated probability of a neutral meson to oscillate first and 
then decay compared to decaying directly is described by the mixing rate 
Rmix- In the absence of CP violation in mixing, Rmix is identical for 
and M^. The following expression then holds for the decay of a into a 

''This convention follows the one used by the CLEO collaboration [7|. 
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Figure 1: The types of processes through which uiixing can proceed in lowest order 
in the SM. Left: Long range processes with ampUtude oc ~iy. Right: Short range 
processes with amplitude oc a: 

final state / [S]: 

i?™, = B{M'> ^W^f)/ B{M'> ^ /) = + y2)/(2 + - (g) 

In the system, the SM expectations for |a;| and \y\ are typically below 
10"'^, which results in RmixiD^) < 10^®. This explains why up to now no 
direct observation of flavor mixing in the system has been achieved and 
why such an observation is likely to be unaccessible also in the future if no 
non-SM effects come into play. 

The situation is much more favorable in the and systems. Experi- 
mentally m, X and —y are known to be of the order unity for the system, 
which implies Rmix{K^) ~ 1. In the system, x ~ 0.7 experimentally PP, 
while y/x ~ —1/500 is expected 0, resulting in R^ixiB^) ~ 0.2.^ 

2.2 Predictions for mixing in the system 

Predictions for x and y both within and beyond the SM span several orders 
of magnitude. For a compilation of predicted values, see Ref. jlj. SM predic- 
tions give typically |x|, \y\ < 10~^ and can be as low as 10~^. Contributions 
of physics beyond the SM may enhance \x\ to values of up to 10^^, while \y\ 
is assumed to be dominated by SM effects. 

At the source of the large uncertainties in the theoretical description, 
even within the SM, is the charm mass being intermediate between heavy 
and light. The charm quark is too light for perturbative treatments to work 

'^For the system, only upper limits have been determined experimentally a; > 
19.0 at 95% CL, while again y/x ^ —1/500 is expected 0, so that Rmix{B^) ~ 1. Direct 
experimental observation of mixing in the system is complicated by the very high 
oscillation frequency oc x. 
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well but, on the other hand, it is too heavy for its decays to be dominated 
by a small number of final states. For a detailed discussion, see Ref. [21 El- 

In the SM, the L>° has the unique feature that its mixing proceeds via 
intermediate states of down-type quarks in the box diagrams of Fig. ^ The 
contribution of the bottom quark to the loops can be neglected because of 
the smallness of jH] Vub x [("i^ — m^)'^ /rny^^m^]. The first factor, Vub, is the 
coupling between quarks of the first and third generation as described by 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix. The second factor corre- 
sponds to GIM-type suppressions in the loop (see below). Consequently, the 

system can be treated to a good approximation as a mere two-generation 
problem that the SM describes with the help of a 2 x 2 rotation matrix. 
There, the rotation angle is the Cabibbo angle 9c, and no CP violating 
parameters are foreseen. 

In that framework, the decay is Cabibbo-favored, while the box di- 
agrams of Fig. fright) are doubly Cabibbo-suppressed. Thus in the 
system, the SM disfavors mixing as expressed by x at least to the level 
sin 0(7^ ~ 0.05. In the SU(3) fiavor limit, where for the quark masses holds: 
niu = rad = rUg, the sum over the box diagrams in Fig. fright) is zero 
(GIM-type cancellation'^) and x would be zero accordingly. In the SM, 
where SU{3) flavor symmetry is broken, GIM-type cancellations result in 
suppression factors that depend on the mass of the heaviest quark in the 
loop relative to the W boson mass: (rrig — m^)^/mp^/m^ ~ rn^/m'^Y for the 
D^, as compared to ~ mf/rriyy (~ Tn'^/m^,) for the {K^)- 

For the on-shell intermediate states of Fig. ^left), which correspond to 
mixing in the system as expressed by y, light quark states are favored by 
phase space. While the limited available phase space for decays results in a 
large lifetime difference in the system (y ~ 1), such constraints for decays 
in the or system are much weaker. The measured branching ratios for 

decays, e.g. into vr'^vr^, K^K'^ or K^n^ are small and not comparable 
to the ones in the system. Based on the same GIM-type suppression 
arguments as for x, SM predictions typically arrive at \y\ < 10~^. However, 
as for the precise size of these GIM-type suppressions is a major source 
of uncertainty. Recent theoretical work |2l E] points out the possibility that 
SU(3) flavor breaking may result in \y\ 0.01 being natural in the SM. Such 
a large value oi\y\ would lower the sensitivity of experimental measurements 
of mixing to potential non-SM physics effects in \x\. 

'^GIM cancellation led Glashow, lUiopoulos and Maiani to infer from the measured 
small size of the ratio BR{K% — > fj,~)/ BR{K1 all modes) the existence of a fourth 
quark, the charm quark. 
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3 Methods for determining x, y in the system 

Given the smallness of the expected mixing rate, Rmix, in the system, the 
extraction of the mixing parameters from a direct observation of flavor os- 
cillations as a function of time is not feasible. Instead, three complementary 
methods have been applied by experiments:'^ 

1. By comparing the lifetimes measured in hadronic decays of specific 
CP symmetry, y can be determined. 

2. Measurements using the wrong-sign decays: 

(a) From the time evolution of wrong-sign hadronic decays, e.g. 
K~^7r~, x'^ and y' can be determined independently, x' , y' are 
related to x, y hy a rotation. 

(b) From wrong-sign semileptonic decays, e.g. — > K^l~i>i, the 
linear combination x"^ + y^ can be determined. 

In order to compare directly the results of method 2(a) to those of meth- 
ods 1, 2(b), an assumption is necessary for the size of the rotation angle (see 
below) between y and x', y' . 

The results of BABAR, Belle and CLEO were obtained with the first 
two methods. They are discussed in detail in the following. Comparisons to 
earlier results by FOCUS (methods 1, 2(a)) and E791 (all three methods) 
are given where appropriate. 

BABAR, Belle and CLEO use for their measurements the following 
hadronic two-prong decays of the D^: 

• The Cabibbo-favored (CF), so-called right-sign (RS) decay — > 

• The CP-even decays D° R- K+ and D° ^ vr 7r+ which are singly 
Cabibbo-suppressed. The former (latter) occurs ~ 1/9 (~ 1/25) times 
less frequently than the RS decay.'' 

• The so-called wrong-sign (WS) decay K~^Ti~ which, in the 
absence of mixing, is doubly Cabibbo-suppressed (DCS) and occurs 
~ 1/300 times less frequently than the RS decay. 

'^The charge conjugate modes are always implied unless stated otherwise. The experi- 
ments treat and D'^ separately only when investigating the possibility of CP violating 
effects. 

''In the limit of SU(3) flavor symmetry, the branching ratios for the decays into 
K~ and tt^tt^ are identical. The experimentally observed substantial difference is an 
example of large SU(3) flavor breaking effects in the D° system [5]. 
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3.1 Measurement of y from a lifetime ratio 

In the absence of CP violation, the final state of the RS decay is an equal 
mixture of CP-even and CP-odd states. Under the assumption that x and 
y are small, the decay-time distribution is approximately exponential [Hj 
with a slope trs = 2/(71 -|- 72), where 71 (72) is the width of the CP- 
even (CP-odd) decay eigenstate Di {D2) defined in Sec. 2. The decay-time 
distributions for the decays into the CP-even final states K^K^ and 
7r~7r"'" are exponential with slope tkk = t-k-k = l/7i- 

Then, in the absence of CP violation, y = Ar/2r = ycp with: 

ycP = 1 or ycp = 1. (7) 

Tkk t't.t, 

In the available measurements y = ycP is assumed.^ 

3.2 Measurement of the time evolution of hadronic WS D° 
decays 

The can arrive at a WS hadronic final state in two ways, either by 
undergoing directly the DCS decay or by first oscillating into a that 
then undergoes a CF decay. This gives rise to three different components in 
the WS decay: from the DCS decay, from mixing and from the interference 
of the decays with (CF) and without (DCS) mixing. 

Assuming CP conservation and expanding the decay rate up to 0{x'^), 
0{y'^) results in the following approximation [3 IHI for the time evolution^ 
of the hadronic*^ WS decay rate: 

Twsit) oc exp(-t) [ Rd + /R^ v' * + 1/4 ( x'^ + y'"" ) t"" ], (8) 

while for the RS decays: Tjisit) oc exp(— t). Here, t is given in units of the 
lifetime of the D^. In this approximation, the time- independent coefficient 
Rd corresponds to the DCS component. Even if the mixing contribution, 
quadratic in t, is very small, the interference term, linear in t, may result in 
a discernible deviation from the exponential time evolution characterizing a 
pure decay. 



'if CP violation is present, then y 7^ ycp and ycp = ycoscp — 0. 5a; A a/ sin 0, where 
Amj^) ^re CP violating parameters related to CP violation in mixing {Am) and in the 
interference of decays with (CF) and without (DCS) mixing (0) [2J. 

^This formula is valid for t <T. 

'^Semileptonic WS Z?" decays can only arise via mixing, so that r^s(i) oc 
exp(-i)( ) t^- 
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The parameters x' and y' are related to the mixing parameters x, y by 
a rotation: 

x' = X cosSk-k +y smdK-K^ y' = y cosSk-k —X sin^xTT- (9) 

The phase, Sk-k, is a strong phase between the DCS contribution and the 
CF one and does not violate CP symmetry. 

In a measurement based on the time-evolution of the WS rate alone, it 
is not possible to determine the phase, Sktz- Furthermore, since only the 
square of x' enters into T\Ys{t)-, the sign of x' cannot be determined, either. 

If there is CP violation in the system, then the parameters Rd, 
x' , y' in Eq.Hhave to be substituted by: R% = ^J{1± Ad)/{1^ Ad)Rd, 
x'^ = \/ (x' cos (piy' sin (p) , y'^ = \/ {y' cos 4>^x' sin <f)) , where = 
(1 ± =F Am)- The plus (minus) sign pertains to the decay of a 

(D^). The three additional parameters are related to CP violation in the 
DCS decay {Au), in the mixing term (Am) and in their interference (phase 

CLEO [7j in its analysis uses a somewhat different approximation which 
is valid for Ad, Am <C 1: P^ = (1 ± Ad)Rd, x'^ = \/l ± Am{x' cos^ ± 
y' sin 0) , y'^ = yj\ it Am {y' cos x' sin (j)) . 

The total time-integrated hadronic WS rate, assuming CP conservation, 

is:^ 

Rws = J Twsit) I J TRsit) = Rd + VR^y' + 1/2 {x'^ + y'^). (10) 

If there is no mixing in the system, then Rws reflects only the rate of 
the DCS decay: Rws = Rd- The SM predicts in the limit of SU(3) flavor 
symmetry: Rd ~ tan^ ^c" ~ 0.0025. SU(3) symmetry breaking effects may 
increase this value (U]. 

3.3 Measurement of the strong phase 6k-k 

No experimental determination of the strong phase, 6ktt, is yet available. 

Reference outlines a method to extract 5ktt from a measurement of 
the rates of the DCS and CF decays of the type D — > Kir. This method 
requires the determination of rate asymmetries in D meson decays to Kltt 
and KsTT- In Ref. Belle finds that the relevant measurements of Kl and 
Ks mesons are possible with a statistical precision sufficient to constrain 
Sr. M- 

'In this approximation, the mixing rate (see Eq.|SJ appears as Rmix = l/2(a;'^ + y'^) ~ 
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4 The experiments 



The experiments BABAR and Belle at the asymmetric B-factories PEP-II 
and KEK-B and CLEO at the storage ring CESR operate near the T(4S') 
resonance at a center-of-mass energy of 10.6 GeV. PEP-II and KEK-B were 
designed with the primary goal of serving as asymmetric B-factories for the 
study of CP violation in the system. In the production cross section, 
a{e'^e~ qq) « 4.45 nb at ^/s ~ 10.6 GeV, bb final states account for jl3j 
~ 1.05 nb and cc final states for ~ 1.30 nb. In terms of production cross 
sections, the B-factories therefore serve equally well as charm- factories. 

Between October 1999 and October 2002, BABAR recorded 94 fb"^ 
of data, while Belle recorded 98 fb~^. The data set analyzed by CLEO, 
recorded between February 1996 and February 1999, corresponds to an inte- 
grated luminosity of 9.0 fb~^. The largest available charm data sets obtained 
in fixed-target experiments were collected at Fermilab by the heavy-flavor 
photoproduction experiment FOCUS and the heavy-flavor hadroproduction 
experiment E791. These data sets contain ~ 120,000 (FOCUS) dH and 
~ 35, 000 (E791) 15 identified — > K~7r~^ candidate events, compared 
to, for example, ~ 260, 000 events n!U| in the BABAR data set. 

Detailed descriptions of the BABAR and Belle detectors can be found 
elsewhere |17( llSj. CLEO carried out its measurement with the CLEO II. V 
detector |19j . Of primary importance for the reconstruction and identifica- 
tion of decays are the vertex and particle-identification detectors. All 
three experiments are equipped with double-sided silicon vertex detectors. 
At BABAR, the Silicon Vertex Tracker (SVT) (five layers, innermost radius 
^min = 3.2 cm) is also capable of stand-alone track reconstruction down 
to particle momenta of ~ 60 MeV. Mounted just outside of the BABAR 
drift chamber (DCH) volume is the ring-imaging detector of internally re- 
flected Cherenkov hght (DIRC). At Belle, the Silicon Vertex Detector (SVD) 
(three layers, r^[^ = 3.0 cm) and the central drift chamber are surrounded 
by an array of aerogel Cherenkov and time-of-flight scintillation counters 
(TOE). The CLEO II. V Silicon Vertex Detector (SVX) had three layers, 
with Tinin = 2.35 cm. 

Differentiating kaons from pions is crucial for distinguishing with high 
purity between the different two-prong decays listed in Sec. 3. CLEO 
relies primarily on kinematic selection cuts based on the measured momenta 
and the assigned mass to distinguish between the different decay modes 
EOl . BABAR and Belle enhance their particle-identification capabilities 
considerably with the help of likelihoods determined from a combination 
of the energy-loss (dE/dx) in the tracking devices and Cherenkov detec- 
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tor information. Based on the ratio of the hkehhoods for the two particle 
hypotheses, BABAR reaches an average identification efficiency (mis- 
identification probabihty) of > 75% (< 8%) for momenta up to 4 GeV pi] . 
Behe includes also the TOF information in the likelihood calculation and 
reports an efficiency of ~ 85% and a mis-identification probability of ~ 10% 
for momenta up to 3.5 GeV |22j . 

5 Measurement of y 

Following Eq. [3 the mixing parameter y is determined by measuring the 
slope of the decay-time distributions in independent candidate samples of 
K~TT^ and the singly- Cabibbo suppressed channels — > K^K^ 
and — > vr^vr'''. 

The selection of a candidate, the determination of its proper decay- 
time and the general structure of the fit to obtain the lifetime from the 
proper decay-time distribution are discussed in some detail in the following. 
They are of importance also for the measurements with WS events described 
in Sec. 6. 

5.1 Method of measurement 

candidates are selected by searching for pairs of tracks with opposite 
charge and combined invariant mass near the expected mass. The com- 
mon vertex of the track pair determines the candidate decay vertex, v^ec-, 
with typical resolutions as listed in Tab. The interception point of the 

momentum vector, p/), with the envelope of the interaction point (IP) 
provides the production vertex, Vprod^ of the L>° candidate [Ml EOl Eg . 

The proper decay-time of a candidate is derived from its mass {mD{PDG) = 
1.864 GeV) ^ and its flight length. In order to take resolution effects prop- 
erly into account, that can, e.g., result in reconstructing a negative value of 
t, the flight length is calculated from the projection of Vdec — Vprod onto po'- 

Ct = mD [ {Vdec - Vprod) ' Pd/\PD?' ]• (H) 

The size of the IP envelope enters into the error on the proper decay-time, 
at, through the uncertainty on the production vertex. The IP envelope is 
smallest in the plane transverse to the beam, as shown in Tab.^ The flight 
length of a candidate is typically ~ 200 ^m in BABAR and Belle. The 
higher decay vertex resolution achieved by CLEO reflects the absence of 
boost-related uncertainties since at CESR, the T(4S') is produced at rest, 
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Table 1: Comparison of parameters that enter into the calculation of the proper 
decay-time, i, at BABAR, Belle and CLEO ^011^1 1 ^ . 





BABAR 


Belle 


CLEO 


IP envelope 








vertical [/im] 


6 


2-4 


10 


horizontal [//m] 


120 


80-120 


300 


along beam-axis [mm] 


8 


3-4 


10 


resolution 


80 ^m 


110 /um rms 


40 ^m 


decay vtx 


along pD 


along p£) 


in each dim. 


resolution of proper decay-time at [fs] 


180 




0.4 X TD 



different from PEP-II and KEK-B, where f3-f = 0.56 and /?7 = 0.425, re- 
spectively, along the beam direction. 

5.2 Event selection 

The three experiments use a similar set of quantities to reduce backgrounds 
jl6L 1201 0^ . The main differences lie in whether a D* tag is required and 
which requirements are imposed for particle identification. 

BABAR and CLEO select events likely to contain a with the help 
of the decay D*^ D^n^, where the vr^ has very low momentum ("slow 
pion", TTg). Each tt^ candidate track is refitted with the constraint to co- 
incide with the candidate production vertex, determined as described 
in Sec. 5.1. This reduces substantially the mismeasurement of the vr^ mo- 
mentum caused by multiple scattering. Then the difference in the invariant 
mass of the D* and the candidate, dm, is required not to deviate from 
the nominal one by more than a certain margin. In BABAR, this margin is 
2 or 3 MeV, depending on whether the tTs reached the DCH; in CLEO the 
margin is 1 MeV. Belle does not use a D* tag and requires solely that the 

candidate flight path be consistent with originating at the IP. 

All three experiments reject events with secondary charm production 
from B meson decays by requiring a minimum center-of-mass momentum 
of the D* or candidates (2.5 GeV for BABAR and Belle, 2.3 GeV for 
CLEO). 

All three experiments reduce background from random combinations of 
two tracks forming a candidate with the help of the angle 6*, measured 
in the candidate center-of-mass system, between the direction of the 

(-D°) candidate boost and its positive (negative) daughter. BABAR 
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Table 2: Data samples used by BABAR, Belle and CLEO in their measurements 

ofy nniisni^. 





BABAR 


Belle 


CLEO 


£ [fb-^] 


57.8 


23.4 


9.0 


K 


vr^ candidate events 


158,000 


214,000 


20,300 




Purity [%] 


99 


87 


91 


R- 


' candidate events 


16,500 


18,300 


2,500 




Purity [%] 


97 


67 


49 




-TT^ candidate events 


8,400 




930 




Purity [%] 


92 




71 



requires cos0* > —0.75 for pion daughters; Belle requires cos 6** > —0.85 for 
pion and, in the K^K^ channel, | cos0*| < 0.9 for kaon daughters; CLEO 
requires |cos^*| < 0.8 in all cases. BABAR and Belle apply in addition 
their particle-identification algorithms to the daughters of the candidate. 
CLEO requires in addition that the invariant mass obtained with the other 
possible particle-type assignments for the two daughters should be more 
than four standard deviations away from the nominal mass. 

The analyses by BABAR and Belle reported here refer only to a sub- 
sample of their full data sets. The sub-sample sizes and the number of events 
after the full event selection are listed in Tab. |2j The purities quoted there 
refer to a window around the mean mass. The window size is ±20 MeV 
in BABAR, ±3cr in Belle {a{K'-K+) = 6.5 MeV, a{K-K+) = 5.4 MeV) and 
±40 MeV in CLEO. Although BABAR's data sample corresponds to twice 
the integrated luminosity of Belle's, Belle's candidate event samples are 
larger than BABAR's and have lower purities because Belle does not require 
a D* tag in its event selection. CLEO arrives at lower efficiencies and, 
for the CP-even channels, lower purities because CLEO, lacking dedicated 
particle-identification detectors, relies primarily on kinematic selection cuts 
to distinguish between the different decay modes. 

5.3 Lifetime fit 

The three experiments determine y from unbinned maximum-likelihood fits 
to the distribution of the reconstructed proper decay-time, t, (see Eq. [TT|) 
of the candidates jl6l 1201 1^ . No background subtraction is performed 
before carrying out the fits; rather, it is left to the fit to describe both signal 
and background. 
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The three experiments proceed in a comparable way to define the full 
fit function: For each of the decay channels, they define a likelihood 
function as the sum of a signal and a background part. Both parts consist 
of a decay-time distribution convolved with a resolution function. The signal 
decay-time distribution is exponential with a slope that corresponds to an 
inverse lifetime. The background decay-time distribution consists of two 
components, an exponential and a delta function, that model background 
events with non-zero and with zero lifetime, respectively. The former can, 
for example, arise from only partially reconstructed three-prong decays; 
the latter corresponds to combinations of particles into a candidate that 
in fact originate at the IP. 

The resolution functions need to accommodate correctly reconstructed 
events and also events with misreconstructed parameters. This need is typi- 
cally met by adding one or more Gaussians to accommodate events that are 
misreconstructed to varying degree. The width of the Gaussian(s) modeling 
the resolution of well reconstructed events is typically given by the per-event 
error on the reconstructed proper time t, at, multiplied with a proportion- 
ality factor. The use of at allows to take into account the uncertainty from 
the highly elliptical shape of the IP envelope (see Sec. 5.1) in the fit. 

All three experiments determine the parameters of the background part 
of the likelihood fit function independently of those in the signal part. Also, 
all three experiments enhance the ability of their fits to differentiate signal 
from background events by multiplying the likelihood function with a per- 
event probability to be a signal event. The further away its reconstructed 
mass, m£), is from the mass peak, the less likely is the candidate to 
belong to the signal. The probability is derived from an independent fit to 
the niD distribution in each decay channel. 

The window in rriD used in the fit must be chosen such that, in addition 
to the signal peak, sufficiently wide sidebands are included so that the 
fit can extract information on shape and size of the background. BABAR 
uses mass windows of 140 MeV width, defined such that reflection peaks, 
e.g., from mis-identified K~tt~^ events in the K~K^ and 7r~7r^ candidate 
samples, are excluded. Belle uses a mass window of ±40 MeV around the 
mean and CLEO uses the window [1.825, 1.905 GeV]. The fit results for 
BABAR and Belle are shown in Fig. |2l The tail at negative values in the 
t distribution reflects the effect of the resolution. Comparing data and fit 
result in the unphysical region t < provides information about how well 
the resolution in the data is modeled by the fit function. The distribution for 
t > reflects the resolution function convolved with the decay-time function. 
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Figure 2: Proper decay-time distribution and fit result. The dots represent the 
data points with statistical error bars. Left: Belle [22]: For Z?" K~tt^ (a) and 
(b), in the D° mass signal region ±3(7 around the mean _D" mass. 
The solid line is the result of the fit, the dashed line indicates the background 
contribution. Right: BABAR jj^l- For D° K~Tr~^, on a linear (upper) and 
logarithmic (lower) scale for all events including the mass sidebands. The solid 
line is the result of the fit, the gray area indicates the contribution attributed to 
background by the fit. 



5.4 Discussion of systematic uncertainties and results 

Belle [5^ determines y from the lifetimes measured in decays into K^K^ 
and K-'K+ {vkk)- BABAR ^ and CLEO ^ use in addition the 'K--K+ 
channel (yTm)- 

Belle and CLEO correct the lifetimes measured in the individual decay 
channels for small biases found with the help of Monte Carlo (MC) studies. 
Belle reports that the proper decay-time tends to be reconstructed slightly 
too small by a decay-channel dependent value. The resulting correction to 
y is (—0.3 lb 0.3)%. The error on this correction, due to the limited MC 
statistics, enters in the final systematic error on the measurement. CLEO 
finds bias values in the individual lifetimes that are compatible, within the 
statistical uncertainty on the MC bias estimate, with zero, chooses to apply 
corrections nonetheless and includes their statistical uncertainty in the final 
systematic error (±1.0% for yxK, ±1-4% for yTm)- BABAR also finds bias 
values compatible with zero and does not apply any corrections to the life- 
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Figure 3: Comparison of results for y from different experiments |14[I15[[T^I2(JII22| . 

The total size of the error bars corresponds to the quadratic sum of statistical and 
systematic error, the inner part indicates the size of the statistical error only. The 
average is calculated as mean of the individual measurements, each weighted by 
the quadratic sum of its statistical and systematic error. 



times but includes in the systematic error on y the statistical uncertainty 
on the MC bias estimate (lo'6% for i/kk, -09% ^or ynTr)- 

Since y is measured from the ratio of lifetimes, many systematic effects 
that affect the individual lifetimes do not affect y. Both Belle and CLEO find 
that the dominant remaining systematic error source is their understanding 
of the background contribution to the signals. In BABAR, the dominating 
source of systematic error is the above quoted uncertainty on the MC bias 
estimate. 

All three experiments arrive at their final systematic error estimate by 
summing in quadrature the contributions of the individual sources. BABAR 
and CLEO each determine y as weighted mean of their measured yxK and 
yiTTT values. 

The results are listed in Fig. 13 Also shown are results from the fixed- 
target experiments^ FOCUS [H] and E791 Within their errors, the 
results are compatible with each other and with the SM expectation of a 
value of \y\ close to zero. 

J The fixed-target experiment results are obtained by means of a fit to the reduced 
proper time. This quantity includes a cut on a minimum distance between primary and 
secondary vertex as a principal tool to reduce non-charm background and is measured 
with a resolution of less than 1/10 of the D" lifetime. The FOCUS (E791) result on y is 
derived from ~ 120, 000 (~ 35, 400) events in the K~n+ channel and ~ 10, 300 (~ 3, 200) 
events in the K~ K'^ channel. 
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6 Measurements with the WS decays K^n , 

DO ^ K-7T+ 

As discussed in Sec. 3.2, the time evolution of the WS decay rate provides 
a means of measuring simultaneously x''^ and y'. CLEO has published a re- 
sult [7] obtained with 9.0 fb^^ of data. BABAR has presented a preliminary 
analysis |5J based on a data sample of 57.1 fb~^. At Belle, work to extract 
the time evolution of WS events is on-going \1'2\ ; results concerning the back- 
ground in the WS sample have already been presented but none yet on 
x''^, y' . All three experiments have results for the time-integrated WS decay 
rate Rws [ZlllIlIlSl- 

6.1 Extraction of x', y' from the time evolution of the WS 
decay rate 

Wrong-sign candidate events of the type — > K^ir^ and D'^ — > K^tt^ 
are selected by requiring the tTs from the D* decay and the daughter K of 
the to have identical charge (WS tag). In addition, CLEO, BABAR and 
Belle apply event selection criteria similar to those described in Sec. 5.2, and 
the proper decay-time is reconstructed as discussed in Sec. 5.1 [7| 1^ OH]. 

The measurement aims at detecting a deviation from a purely exponen- 
tial decay law in WS events (see Sec. 3.2) by means of a likelihood fit to 
the distribution of the reconstructed proper decay-time, t. The likelihood 
functions used by the experiments have the same principal structure as the 
ones used in the measurement of y (see Sec. 5.3). They differentiate between 
a signal and a background component and model each as the convolution 
of a decay-time distribution and a resolution function. For WS events, the 
decay-time distribution follows Eq. |H1 i.e. includes the three parts discussed 
earlier: DCS decay, mixing and interference between the decays with (CF) 
and without (DCS) mixing. 

The WS sample has only ~ 1/300 times the number of events of a RS 
sample selected with the same criteria as the WS sample except for the WS 
tag. This RS sample is used wherever possible to constrain aspects in the 
fit that are common to the two samples. CLEO and BABAR determine, 
e.g., the resolution functions of WS and RS signals and of the common 
background types with the RS sample. 

In the WS sample, a significant additional complication arises from the 
much lower achievable purity'^, e.g. ~ 50% in the CLEO analysis. The fit 

''A rough estimate can be obtained by assuming that in a given random sample of 
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D*+_^ D °7r+, D°-^K'^7r" 3280ioo-ooi 




Figure 4: Signal (dots, with statistical error bars) observed by CLEO 7 for the 
WS decay _D" K^n^ (charge conjugate mode included). The projections of 
the fit for the signal (cross-hatched) and the background (background of type A 
- "uds", type B - "charm", type C - if+Tr-") is shown. M (left) and Q 

(right) are each within 2a of their mean in the RS mode. 

needs to accord each background type its specific lifetime evolution, distinct 
from the ones of the other background types and from those of the signal 
events. The reconstructed mo and 6m distributions provide a means of 
differentiating between the different background types and between back- 
ground and signal. 

6.1.1 Fitting procedure 

CLEO performs a two-step fitting procedure. First, the levels of the dif- 
ferent background types in the selected candidate sample are estimated by 
performing a fit to the two-dimensional region 1.76 < uid < 1-97 GeV and 
< Q < 10 MeV, where Q = 5m — mj^ and is the pion mass. CLEO dif- 
ferentiates three types of backgrounds: A) from the CF decay that together 
with a random vr^ mimics a WS event ("Z)'^ K~^tt~"), B) with non- 
zero lifetime ("charm") and C) with zero lifetime ("uds"). The background 
shapes in the — 5m plane are estimated with the help of a Monte Carlo 
sample that corresponds to £ = 90 fb~^. The signal shape is constrained 
with the help of the RS sample. The fit yields ~ 45 WS signal events in the 

events, the probability for a non-signal event to be mis-identified as a signal event is 
the same for WS and RS signal. Then, a purity of 99% in the RS candidate sample 
translates into a purity of 25% in the WS candidate sample because WS events are 300 
times less frequent than RS events. This simplified example does not include background 
from random ttb that together with a CF decay mimic a WS event. 
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Figure 5: (a) Distribution of the proper decay-time, t, in the WS data sample 
(dots) compared to the fit result for the different background types (A: light grey, 
B: black, C: darker grey) and for the signal (open), (b) 95% C.L. contours, including 



the systematic uncertainty, separately for (y'+ vs. a;'"*"^) and {y' 



vs. X 



when allowing for CP violation, and 13°, combined {y' vs. x'^) when not allowing 
for CP violation. The solid (open) point represents the most likely fit point in the 
case of no CP violation without (with) the constraint x'"^ > O.pT) 



total of 82 selected candidate events, measured within a 2a window around 
the central values of the rriD and Q distributions, (a (mo) = 6.4 it 0.1 MeV 
and cr{Q) = 190 it 2 keV). Figure E] shows the observed signal as function 
of Q and M = uid together with the projected background estimates from 
the fit. In the second step, CLEO performs a binned maximum-hkelihood 
fit (bin size 1/20 of the lifetime) to the distribution of the proper decay- 
time, t. In this fit, the background levels from the first fit in the tjid — Sm 
plane are used. 

BABAR performs an unbinned maximum-likelihood fit in several steps. 
First, the number of signal and background events is determined in a fit 
to the m£) — 5m plane. In the WS sample, BABAR models separately the 
background from A) the CF decay that together with a random vr^ mimics 
a WS event, B) RS events where the K and vr hypotheses are swapped, C) 
purely combinatorial background. In the RS sample, BABAR considers in 
addition background from D) partially reconstructed mesons, typically 
from 3-prong decays. The 6m distributions for background types A) 
and C) are obtained directly from the data by means of mixing a vr^ from 
a D* decay of one event with the candidate of another event. The fit 
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yields ~ 440 signal WS events in the region 1.804 < rriD < 1.924 GeV, 
6m < rriTr + 25 MeV. Figure IHIa) shows the reconstructed proper decay-time 
distribution for the WS signal and the three background types A)-C). In 
a second step, a fit to the RS and WS proper decay-time distributions is 
performed simultaneously to determine the resolution functions. 

Belle has reported results on the first part of the fitting procedure, the 
two-dimensional fit in the rriD — 6m plane. Belle considers background types 
similar to BABAR's. Specifically to suppress background type B), Belle 
requires that the invariant mass of the candidate calculated with swapped 
particle hypotheses for the daughters differs from the nominal mass by more 
than 28 MeV (~ 4cr). The remaining background types are modeled in the 
mD — 6m plane fit function. The fit arrives at ~ 450 WS signal events in the 
region 1.81 < mo < 1-91 GeV, < Q < 20 MeV. Belle has not yet released 
any results from the fit to the distribution of the proper decay-time, t. 

For the WS decay-time function, CLEO and BABAR both follow Eq. H 
and consider the case with and without CP violation. Input to all fits is 
mD, 6m, t, at. 

6.1.2 Discussion of results 

The results of both CLEO [7j and BABAR ^ are consistent with the 
absence of mixing and of CP violation. 

CLEO constructs contours in the x' — y' plane that contain the true value 
of x', y' with 95% confidence. Following a Bayesian P approach, the 95% 
confidence region includes all points in the plane for which the negative log 
likelihood, — InL, differs from the best fit value by more than 3.0. All fit 
variables other than x', y' are allowed to vary to give the best fit value at 
each point on the contour. The 95% confidence region, when not and when 
allowing for CP violation, is shown in Fig. Efa). There, the assumption 
8 Kit = 0° is made for the strong phase (see Eq.lHJ, so that x' = x and y' = y. 
The contours do not include any estimate of systematic errors. CLEO also 
quotes one-dimensional limits at the 95% C.L. for the mixing parameters, 
(l/2)x'2 < 0.041% and -5.8% < y' < 1.0%, without assumptions on CP 
violation. CLEO identifies as its dominant source of systematic error the 
potential misunderstanding of the background shapes and acceptances and 
assigns a total systematic error of ±0.2% (±0.3%) for x' {y'). 

Figure Efa) also shows a preliminary contour from FOCUS |2n|, derived 
from an analysis of WS hadronic decays of the type — > JC+vr", a 95% 
C.L. contour from E791 derived from the 90% C.L. limit Rmix < 0.50%^ 

'Rmi. ~ l/2(x'^ + t/'2) = l/2(x^ + y^). 
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Figure 6: (a) Comparison of the 95% C.L. contours in the x — y plane from 
E791 ,2^, FOCUS 26^ and CLEO 0. For the bigger kidney-shaped (smaller) 
region from CLEO, CP violation was (was not) allowed in the fit. Also shown is 
the band corresponding to the current average for y obtained in lifetime-difference 
analyses (see Fig. O. The assumption 6ktt — 0° for the strong phase is made, (b) 
Comparison of the 95% C.L. contours in the x'^ — y' plane from CLEO and BABAR 
when assuming CP conservation in the fit. The CLEO contour docs not include 
systematic uncertainties. 



in WS semileptonic decays of the type D K^l^vi, and the band corre- 
sponding to the average on y from the lifetime difference analyses (see Fig.E)) . 
Not shown is the contour from E791 [21] derived from the weaker 90% C.L. 
limit Rrnix < 0.85% in WS hadronic decays of the type —>■ K^Tr~ and 

BABAR uses a method different from CLEO's to arrive at a 95% C.L. 
contour. BABAR argues that due to allowing x'^ in its fit to take unphysi- 
cal negative values, it is not clear how to apply a Bayesian ansatz to derive 
an error estimate from the two-dimensional likelihood distribution. In ad- 
dition, BABAR finds that the likelihood distribution depends strongly on 
the most likely fitted values of x'^ and y' . To define 95% C.L. contours, 
BABAR applies a frequentist approach based on toy Monte Carlo (MC) ex- 
periments. Any point a*c = {x'c^y'c) on the 95% C.L. contour has to meet 
the requirement: If a toy MC experiment is generated at that point, there is 
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Table 3: Comparison of the 95% C.L. limits for the fit output parameters of 
BABAR and CLEO when CP conservation is assumed in the fit. BABAR's limits 
include systematic uncertainties and were obtained in a fit that allowed x'"^ < 0. 





Rd 


y' 


x' 


x'^ 


CLEO [%] 
BABAR [%] 


(0.24, 0.69) 
(0.22, 0.46) 


(-5.2, 0.2) 
(-3.7, 2.4) 


(-2.8, 2.8) 


< 0.076 

< 0.21 



a 95% probability that the ratio A In L(d*c) = InL(crc) — In Lmax is greater 
than A In Ldatai^'c) calculated for the data. There, L^ax is the maximum 
likelihood obtained in the fit either to data or to a toy MC sample. BABAR 
constructs 95% C.L. contours also for the systematic effects considered in the 
analysis, among them uncertainties in the form of the fit functions, detector 
effects and effects of the selection criteria. Figure E^b) shows the resulting 
95% C.L. contours that include the statistical as well as the systematic un- 
certainty estimate. A strong correlation between x'"^ and y' is apparent. The 
most likely fit point in the case of no CP violation has a negative coordinate 
m X . 

A direct comparison of the CLEO and BABAR results is not possible for 
the case when CP violation is allowed in the fit. CLEO uses as fit output 
parameters x' , y', Rd and Ad, Am, sincp, while BABAR uses x'"*"^, R~^ 
(x'-2, y'-, R-) for the (W) case (see Sec. 3.2). 

The results when assuming CP conservation in the fit are in principle 
comparable between CLEO and BABAR, see Tab. 01 BABAR includes 
systematic uncertainties in its 95% C.L. contour and obtains its limits on 
x'"^ and y' by projecting this contour onto the corresponding axis. The CLEO 
limits in Tab.|31 instead, correspond to one-dimensional 95% C.L. intervals, 
determined by an increase in — InL of 1.92 compared to the best fit value, 
and do not include systematic uncertainties. 

BABAR's upper limit on x'^ is almost three times bigger than CLEO's, 
in spite of being based on a six times larger data sample. A possibly overly 
conservative estimate of the systematic error cannot account for a difference 
of this size, as illustrated by Fig. E^b). There, the CLEO 95% C.L. contour 
(no systematic errors) is overlaid with the two BABAR contours that are 
obtained before and after adding the systematic uncertainty. Two possible 
reasons for so pronounced a difference between CLEO and BABAR are the 
different techniques for obtaining the 95% C.L. limits and the treatment of 
the fit output parameter x'. BABAR allows x'^ to take unphysical negative 
values in its fit, while CLEO excludes this possibility by choosing x' instead 
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Figure 7: Comparison of results for Rws from different experiments [3 1211 
1241 1281 ■ The total size of the error bars corresponds to the quadratic sum of 
statistical and systematic error, the inner part indicates the size of the statistical 
error only. The average is calculated as mean of the individual measurements, each 
weighted by the quadratic sum of its statistical and systematic error. 



of x'^ as fit parameter. Another possible reason is the sign of the fit result 
on y' . Toy Monte Carlo studies indicate that the 95% C.L. contour differs in 
size and shape depending on the sign [I'll I27j . For positive y' (the BABAR 
case), the contour tends to be larger than for negative y' (the CLEO case). 
Albeit its precise origin is not yet understood, this behavior may point to a 
qualitative difference between the two regimes, with destructive (construc- 
tive) interference between the decays with (CF) and without (DCS) mixing 
for y < (y > 0). 

6.2 Extraction of the time-integrated WS decay rate Rws 

CLEO and BABAR arrive at a measurement of Rws by repeating the fits 
described above with the assumption of no mixing in the system, i.e. 
x = y = [Tlin]- 

Belle uses its fit in the mD — Sm plane, described above, to determine the 
time-integrated number of signal events in the candidate samples |23| . The 
ratio of the number of signal events in the WS and RS candidate samples 
yields Rws- The systematic error on Rws is dominated by the uncertainty 
on the background shapes used in the fit. 

Figure [71 compares the results from CLEO [7j, based on £ = 9.0 fb~^, 
BABAR (21] (57.1 fb"^) and Belle |2H] (46.2 fb"^) to earlier measurements 
by E791 HU, ALEPH 28_ and FOCUS ^ in WS decays of the type D° ^ 

K+TT-. 



24 



In the absence of mixing, Rws corresponds to the the doubly Cabibbo- 
suppressed decay rate Rd (see Sec. 3.2). The discrepancy between the SM 
expectation in the absence of mixing, R^) 0.25%, and the values of Rws 
measured by CLEO, BABAR and Belle may be attributable to the effect of 
SU(3) symmetry breaking [S]. 

7 Summary and Outlook 

The asymmetric B-factories PEP-II and KEK-B with the experiments BABAR 
and Belle, operational since 1999, have rendered possible studies of mixing 
in the system with unprecedented statistical precision and sample puri- 
ties. First results for the mixing parameters x and y are compatible with 
an absence of mixing and of CP violation. 

BABAR and Belle both expect to reach C = 500 fb"^ by 2006. In 
data sets of this size, the statistical errors in the measurements reviewed 
in this article can be reduced by a factor of three. For y, a statistical 
precision of ~ 0.2% would be within reach and could be further improved by, 
e.g., employing additional decay channels of well-defined CP symmetry 
and, for BABAR, by investigating the possibility of dropping the D* tag 
requirement in its event selection. 

BABAR's preliminary result for x'^, y', determined from the time evo- 
lution of WS decays of the type K~^7r~, is the second such result 
from a collider experiment after CLEO's. Already with its present data 
sample, BABAR should be able to improve its systematic error estimate 
substantially. Future results from Belle should help shedding light on the 
differences between the CLEO and the BABAR result. Additional infor- 
mation may be gained from other WS decay channels, hadronic and 
semileptonic ones, which have already been observed in other experiments 
(see e.g. Ref. [211120]). 

Given the huge uncertainties in Standard Model (SM) predictions for 

mixing, it might prove difficult to establish physics beyond the SM from 
the size of the measured x and y parameters alone. A more robust potential 
signal for new physics may well be CP violation in the system ^ . Efforts 
to establish CP violation in the system are already part of the WS time 
evolution analyses discussed in this review and are likely to intensify in the 
future. 

Methods for investigating mixing in the system complementary to 
those at the B-factories would be available at CLEO-c. CLEO-c at CESR-c, 
currently under discussion as successor to CESR at Cornell, may operate 
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at the D^D^ threshold. There, the quantum-mechanical coherence of the 
produced D^, pair can be exploited to study mixing in the system 
in ways not available™ to any already existing experiment [311 132j . 
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